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Abstract 

A finite-size correction to a probability to detect the neutrino pro- 
duced in high-energy hadron reactions is studied in this series of papers 
I and II. The neutrino is produced as a secondary particle from a weak 
decay of stable particles produced in hadron reactions and detected 
at a distant position with a target nucleus. In I, we describe a general 
formalism on the finite-size corrections and particularly on those of 
the pion in hadron collisions. The finite-size correction has an origin 
in a many-body wave function at a finite time in which the kinetic 
energy is not a good quantum number. The state is non-uniform in 
space-time and reveals interference pattern that can be computed rig- 
orously via the light-cone singularity of the involving particles. The 
correction depends on the mass and energy of the observed particle 
and becomes small for the pion. In the following paper II, the neutrino 
in the pion decay is shown to reveal a large finite-size correction that 
depends on the absolute neutrino mass in a universal manner. Thus 
a new method of measuring the absolute neutrino mass is suggested. 
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1 Introduction 



Understanding neutrinos is in rapid progress these days. Neutrinos are very 
light and mass squared differences were given from recent flavor oscillation 
experiments [U El El HI EJ E] using neutrinos from the sun, accelerators, re- 
actors, and atmosphere as [7], 

Am 2 21 =m\-m\ = (7.50 ± 0.20) x 10" 5 [eV 2 /c 4 ], (1) 
\Amj 2 \ = \ml - m\\ = (2.32±°3) x 10~ 3 [eV 2 /c 4 ], (2) 

where m» (i = 1 — 3) are mass values. Absolute masses are important but 
are not found from oscillation experiments. Tritium beta decays [8] have 
been used for determining the absolute value but an existing upper bound 
for an effective electron neutrino mass-squared is of the order of 2 [eV 2 /c 4 ] 
and the mass is 0.3 — 1.3 [eV/c 2 ] from cosmo logical observations [9]. In 
these neutrino experiments, higher precision and more statistics have been 
achieved and will be improved more. In this series of papers, we study an 
interference phenomenon of neutrinos that could give informations of the 
absolute masses. 

A neutrino interacts extremely weakly with matter so is idealistic for a 
test of quantum mechanics. Although its validity is certain from many tests 
and applications made in the past, most of them are restricted mainly to 
microscopic areas. Neutrinos provide a new test in macroscopic regions. For 
the test, the number of neutrino events must be large to have small sta- 
tistical fluctuations. In electron bi-prism experiments by Tonomura et al 
[TO] , interferences of electrons of a sharp energy become visible when the 
total number becomes significant. Even though initial electrons are created 
randomly, when the number of events becomes large and a signal exceeds a 
statistical fluctuation, a distribution of detected electrons displays an inter- 
ference pattern. Each pion decays randomly but the entire decay process is 
described by a quantum mechanical wave function and the distribution of 
the produced neutrino could display a similar interference pattern. For the 
interference pattern to be observed, a large number of events is necessary. 

Charged current weak decays occur with an interaction Lagrangian of a 
product of V — A currents in low energy region and a neutrino and charged 
lepton are produced in pair. Weak decay processes have been studied using 
an S-matrix of plane waves with asymptotic boundary conditions, where 
particles in the initial and final states are regarded as free waves without 
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correlations [TTJ [TJl EH [III EE]- A decay rate, average life time, and various 
distributions of charged leptons have been computed, and perfect agreements 
with experiments have been obtained [T6]. The charged leptons are in fact 
in the asymptotic space-time regions in these experiments. Probabilities to 
detect particles in asymptotic regions do not depend on distances between 
positions of the production and the detection, and are computed with the 
standard S-matrix of plane waves. This method may not work for neutrinos, 
since neutrinos propagate with almost the light speed due to their small 
masses and are detected at a distant position. Hence the neutrino wave 
produced in a pion decay is similar to an electrostatic potential of a moving 
body, which has a finite-size correction in a form of a retarded potential. 

A wave function that starts from an eigenstate of a free Hamiltonian H 
of an eigenvalue Eq and evolves with a total Hamiltonian Hq + H\, where 
Hi is an interaction Hamiltonian, is a superposition of waves of various ki- 
netic energies Ep at a finite time. So the wave function expressing the whole 
process approaches an asymptotic form slowly in certain situations. For com- 
puting a probability in this region, the ordinary method of S-matrix, which 
we write S'foo], is useless and those that satisfy the boundary condition at a 
finite-time interval T and wave packets localized in space that fulfill bound- 
ary conditions of experiments are necessary. An S-matrix S[T] that satisfies 
the boundary conditions at a finite T is appropriate to study them. We de- 
velop a theory for S[T] and compute the finite-size corrections of transition 
probabilities with S[T] expressed by wave packets. 

An amplitude to detect a particle at a finite distance or time is written 
with suitable wave packets for the final states in LSZ [14J reduction formula. 
When calculations are made without replacing the wave packets with plane 
waves, an expression is obtained and the probability is computed. To study 
the probability that depends on the position is unfamiliar in high-energy 
physics, hence we review an essence of the method first. Let ip{x, 7> X^) 
be a wave function at a space-time position x of a particle that we observe 
and other particles 7 in a scattering process, where X® shows the position 
of the target. ip(x, 7; X^) is expressed by a corresponding Green's function 
and its explicit form for a detection of a pion will be given later and that 
for a neutrino will be given in II. The amplitude to detect the particle of a 
momentum p at X is written as 

*(X,p, 7 ) = / d 4 xw(p,x-X)^x, Tl X^), (3) 
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where w(p,x — X) is a wave packet that satisfies a wave equation of the 
particle and X® is omitted in ^(X,p, 7). A Gaussian form with a size a of 
a momentum pat a position X^ 

W{p, X-X)= Noe -(S-X-*(t-X°)r/<r-i P -(x-X) } (4) 

is used for the sake of simplicity, where Nq is the normalization factor and v 
is a velocity of the relativistic particle Vi = -J^-E(p) = p { /E(p). We write a 
four- dimensional coordinate as x and its components as (x°, x). For the 0-th 
component a; , sometimes t is used, x° = t. 

An integral over the states 7 of the square of modulus of this amplitude 

C(X,p) = J d7**(X,p;7)*(X,p; 7 ), (5) 

is a probability of detecting this particle at X, where four dimensional coor- 
dinates x are integrated over < x 11 . X® is the coordinate of the target 
particle, and expresses the whole phase space of the state 7. 

C(X,p) is written with a correlation function A (07, £2) defined from the 
wave function ^(£,7; X^) as 

A(x u x 2 ) = j d T v>*(zi,7;^ (i) )V^2,7;* ( ^ (6) 

This function, in a relativistic invariant system, has a term of the form, 

d P J(p+Q)-(xi-x 2 ) ( 7 \ 
E(pf ' [() 

where p M is the four- dimensional momentum of the mass shell of a particle 
and is a momentum of others. This integral is known to be composed of 
a term proportional to the singular function [T7| 

6(X)e(t 1 -t 2 ),X = (x 1 -x 2 ) 2 , (8) 
e{h - 1 2 ) = ±1 for h-t 2 ^ 0, (9) 

and those proportional to the regular functions e* mA , e~ mX where a finite m 
is determined from a dynamics of the system. 
C(X,p) is expressed with A( X\ 5 X2) 8jS 



C(X,p) = / Y[ d 4 XiW*(p,xi - X)w(p,x 2 - X)A(xx,x 2 ), 

J i=l,2 



(10) 
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and behaves as 



C® + CWg(\X-X®\/lo), (11) 

where and depend on the momenta and a function g(\X — X^\/l Q ) 

is studied in Section 3 and decreases rapidly in ^ X ~^ - > 10, and behaves 

as -g 2 % ru at \ x ~^ { } l — >. oo. The rapidly oscillating or decreasing terms 

in A(xi,X2) contribute to the constant term, C^°\ and the singular term, 
Eq. ([!]), gives the second term of the above equation expressing a finite- 
size correction. l determines a typical length that the finite-size correction 
remains and is called a coherence length of the wave function. Thus the 
singular function gives the long-range correlation. 

From explicit calculations which are presented later, a form of Zo m a high 
momentum region, \p\ 3> m, is given by 

where m is the particle's mass. For a pion and electron of energy 1 [GeV], 
they are 

C° n = ^[GeV- 1 ] = 2xlO- 14 [m], (13) 

electron = jf^^y-l] = ± 3 x 1Q -9[ m ]. ( 14 ) 

Other hadrons are heavier and have shorter lengths than that of the pion. 
Thus Zo for hadrons and charged leptons are microscopic lengths. 

In the asymptotic regions where initial and final states are expressed 
by free particles, C(X,p) neither depends on the longitudinal coordinate 
(X — X^')l parallel to the momentum p nor the time X° - X o(l) and the 
wave behaves like a particle and the standard S-matrix is applied in this 
space-time region. Processes of detecting particles in \X — X®\ Zo are in 
the asymptotic region. Since the particle in this region is expressed by a plane 
wave and behaves like a free particle, the particle number and flux behave 
like those of classical particles and the production and detection are treated 
separately. The particle flux is determined by its distribution functions of the 
number and velocity that are determined by the decay process and is used for 
calculation of the scattering processes. From Eq. (|14p . scattering experiments 
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Fig. 1: The whole process in a high-energy neutrino experiment is illustrated. 
By a collision of a proton with a target, a pion is produced. The pion 
propagates a macroscopic distance and decays in a decay tunnel. A neutrino 
is produced and is detected. 



of the charged leptons and hadrons using detectors of macroscopic sizes and 
distances are treated with the ordinary S-matrix. The transition amplitude 
has a delta function of an energy and momentum conservation and a total 
probability becomes proportional to a large time interval T and a decay 
rate and a mean life time are found from the probability per time, i.e., a 
proportional constant. 

If C(X,p) varies with (X - X®) L or X° - X° [i) , the final state behaves 
like a correlated wave and is not in the asymptotic region. The process then is 
not described by the standard S-matrix of plane waves but by a Schrodinger 
equation of a suitable boundary condition, which describes a time-evolution 
of the state vector regardless of the boundary conditions, or by the S-matrix 
expressed by wave packets and satisfies a boundary condition at a finite T, 
S[T}. 

The whole process of the production and decay of a pion is expressed 
in Fig.[TJ The first reaction is caused by strong interactions and the second 
reaction is caused by weak interactions. These interactions are understood 
well and scattering probabilities, cross sections, and other physical quantities 
which are defined at asymptotic regions and depend on particles' momenta of 
the initial and final states have been computed theoretically and agreements 
with experiments have been obtained. 

Here we study new possibilities that the particles are detected in non- 
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asymptotic regions where C(X,p) varies with (X — X^)l or X° — X ^. 
The pion is studied in I and the neutrino is studied in II. We find that 
the pion retains wave nature in short distance region and is detected in 
the asymptotic region but the neutrino retains wave nature in a wide area 
and is detected in the non-asymptotic space-time regions in certain area. 
A probability to detect a neutrino varies with the distance and has a large 
finite-size correction, in fact. Being described by a many-body wave function 
together with the charged lepton and pion, the neutrino is not a simple 
isolated wave but is a superposed wave of many components at a position 
of detector. The neutrino is detected through its interaction with a nucleus 
in a target. Hence the probability to detect neutrinos at a finite distance is 
computed with wave packets that are localized in space. The wave packets 
vanish at a position x if the distance between x and the center position X, 
\x — X\ — > oo, and satisfy the boundary conditions of the experiments. Hence 
they are appropriate to study the finite-size effect and are used here [HJ [T5] . 

This paper is organized in the following manner. In section 2, S-matrix 
of a finite-time interval, S[T], is introduced. In section 3, pions in hadron 
reactions are studied. Summary and prospects are given in section 4. Various 
properties of wave packets including the size, shape, and completeness are 
studied in Appendix. The neutrino in the pion decay is studied in II. 

2 Finite-size correction and S-matrix at a finite- 
time interval 

In a scattering amplitude, an initial state is prepared at t = — oo and a final 
state is observed at t = +oo. In practical situations, these times are not 
infinite but large values ±T/2. Approximation of the finite T with oo is 
normally good and a deviation of physical quantity between those of T and 
of oo, which we call a finite-size correction, is negligibly small then. That is 
not the case in some situations. So, general properties of wave function and 
scattering matrix at the finite-time interval are studied, first. 



7 



2.1 Wave function at a finite time 
2.1.1 Wave functions at a finite time 

In a physical system described by a Hamiltonian H composed of a free term 
H and an interaction term H 1 

H = H + H l7 (15) 

the Schrodinger equation in a unit of h = 1, 

^|*(t)) = (ffo + ffi)|*(t)) ( (16) 

describes the time evolution. Unitary operators 

U(t) = e~ iHt , U (t) = e~ iHot , (17) 

give time evolutions of the state vectors. State vector in the interaction 
representation defined by 



satisfies 



|*(*)) = U \*(t)), H^t) = UotyH&iit), (18) 



and a solution is given by a time-ordered product, 

\V(t)) = T f dt'e^'^iO)) = |tf(0)) + f dt'A(t')\V(0)), (20) 
Jo Jo 

where 

A(t') = H 1 (t')/i+ f dt"(H 1 (t')/i)(H 1 (t")/i + ---). (21) 
Jo 

In the first order in Hi, there is no ambiguity of higher order corrections 
and the solution is 

|¥(*)> = \l + j^dt'Hi{t')/i} |¥(0)>. (22) 



and is written as 



|*(0))+ I dt' I dp\P){P\H x {t')/i\*{<S)) 



/ lU)t _ 1 
dP |0></9|#i(O)|*(O)), (23) 
u 

where 

#o|* (0)> = Eo\*(0), H \P) = E \/3), (24) 

and 

u = E p - E . (25) 

At t — > oo, the formula 

e iuJt - 1 fwx(ud/2)' 



2ie^ — ±-!—L ~ 27ri(5(w), (26) 

UJ \ CO J 

is substituted, and the states of Ep = E contribute. But at a finite, t 
states of broad spectrum contribute. At a finite u, the exponential function 
oscillates rapidly and its average in St of satisfying u5t ^> 1 is given by 

G - 1 » -i. (27) 

a; a; 

The wave function Eq. (12 3 p at the infinite t and the average at a finite t agree 
with 



^(^oo) = |*(0)) + 27ri / dp\l3)(P\Hi(Q)WW(Ep-E ), (28) 

|*(*)o«en*e) = |#(0)> - / dfi- \0) (/3|#l(0) | #(0)> J finite t. (29) 

Both of them have the total energy E of the total Hamiltonian H 

H\V(t)) = E \y(t)). (30) 
They have different energies of Hq 

H \^(t) oo ) = E \^(t) oo ) 1 (31) 

Hoy^{tj average) 7^ ^0\^ (t) average) ■ (32) 
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It is noted that the energy Ep of H in Eq. (|2"5|) is E and the kinetic 
energy defined by the free part H is conserved at the infinite time, whereas 
Ep in Eq. (1231) is in a wide range of certain weight and a mean value of Ep 
varies with time and differs from Eq. The difference is due to H\. Hence the 
energy defined by the free part Hq is not conserved at a finite time but the 
difference Ep — E vanishes in the asymptotic region. 

Thus states at t — > oo are governed by the states of the same kinetic 
energy Eq and the asymptotic values such as the total cross section and 
decay rate are computed within this space. Now the state at a finite t, 
Eq. (1231) . is completely different and is composed of many components of 
eigenvalues Ep ^ E . This state is a superposition of waves of different 
kinetic energies with a definite weight and is similar to electron waves formed 
by interferometer such as that in a bi-prism experiment. Plane waves of 
various kinetic energies are generated by potential energies of an obstacle or 
potential in a stationary wave and show a clear single-electron interference 
phenomenon. 

Here waves of various kinetic energies are generated by the interaction 
energy due to Hi in a time-dependent wave function. This many-body in- 
teraction initially makes a transition of a particle to two particles if that is 
a form of 3 or from two particle to two particle if that is a form of 4 and 
gives an energy to a state at a finite t of Eq. (f23|) as a by-product. Since the 
total energy is conserved, the kinetic energy is not constant. Accordingly, 
the wave Eq. (1231) could show a diffraction phenomenon that is characteristic 
of the waves of many components of certain weight. This diffraction pattern 
depends on t and gives a finite-size correction to the probability. They de- 
pend upon the spectrum and states at all Ep, even this is the phenomenon 
of the wave function at the tree level. 

The weight of \/3) in wave functions at a finite t and t = oo are pro- 
portional to e "" ~ 1 and 2iri5(cu) in any order of Hi. Thus the waves are 
superpositions of waves of all Ep at a finite t, and superpositions of waves of 
Ep = E at t = oo. 

2.1.2 Scattering operator of a finite-time interval 

The wave function is not directly observable. Observables are measured 
with scattering or decay processes of the particles and their amplitudes are 
described by M0ller operators. Those at a finite-time interval T are defined 
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from U(t) and Uo(t) in the form 

Ofc(T) = lim lP(t)U (t), (33) 

t— >=fT/2 

and satisfy 

e iHr fi T (T) = fi T (T ± r)e^° r , (34) 
Taking an infinitesimal r, we have 



1 + i#r)Q T (T) = |o T (T) ± r^O^T)} (1 + M t). 



(35) 



Thus the relation 



HQ T (T) = T i Aq t ( T ) + fi T (T)# , (36) 

is satisfied. The first term in the right hand side is a specific term of those 
operators of the finite T. Hermitian conjugate of this equality is 

Q^(T)H = ±i (±t±(T)j + H°0^(T). (37) 

Scattering operator of a finite-time interval T is the product 

S(T) = f)L(T)fl + (T), (38) 

and satisfies 

S(T)H = nl(T)Q + (T)H 

= fi t _(T)|iJfi + (T)-^ fi+ (T)} 

= jy nL(T) + i (^?n-(T)) n+(T) -^L(t)Aq + (t) 



# s(t) + < ^n_(T)j n+fT) -inLcrjAn+cD, 



(39) 



from Eq. ( 136 p . Finally we have 

[S(T) , H ] = i (^«- (T)) fi + (T) - iOL (T) (T). (40) 
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S(T) does not commute with H and does not conserve the energy defined 
by H . The matrix element has the conserving and non-conserving terms, 

(P\S(T)\a) = 8p, a + 5 e (E a - Ep)f a> p + 5f(T) (41) 

of the kinetic energy. Since the energy Ep of the second term of the right-hand 
side of Eq. (14ip is different from that of the third term, the total transition 
probability becomes a sum of each probability. The second term gives con- 
stant physical quantities that are computable by the ordinary S-matrix of 
plane waves and the third term gives T-dependent corrections that are not 
computable by the ordinary S-matrix. A magnitude of Sf depends on a dy- 
namics of the system and a probability gets modified by this term at the 
finite T. When E a and Ep are approximate energies of the states |a) and 
we have 

(E a -E p )(P\S(T)\a) = (P\0(T)\a), 

o(t) = i (^n-(T)) n + (T) - <nL An+cr). (42) 



Hence 



Sf=- -(0\O{T)\a), (43) 



and a probability from these states that have the energies Ep different from 
E a is given in the form, 

2 |(/3|O(T)|a)| 2 >0, (44) 



where the equality is satisfied at T — > oo. Thus the finite-size correction is 
derived from those states that have the energy Ep ^ E . The examples of 
the finite-size correction in the lowest order of Hi that reflects the states at 
ultra-violet region will be presented later in I and in II. 

The fact that the energy defined by a free part H Q is not conserved in a 
transition process of a time interval T, even though the total energy defined 
by a total Hamiltonian H is conserved and the state /3 of the energy Ep 
satisfying Ep 7^ E a contributes is known well in quantum mechanics [18]. A 
probability of detecting a state /3 at T is given in the form, 



2 4sm a [{E p -E a )T\ 



(Ep — E a 
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where F a> p is the matrix element. A variance of the energy difference AE = 
E a — Ep satisfies an uncertainty relation 

TAE > h. (46) 

In a system of continuous spectrum also, a formula 

is used normally. If a system has states of \Ep — E a \ = e, then the uncertainty 
relation leads the time T must satisfy 

T > 5. (48) 

Conversely at a time below this value, the amplitude and probability may 
be different from the asymptotic values, and may have large finite-size cor- 
rections. The large T behavior of Eq. ( I45j) should be examined carefully and 
not be replaced with the asymptotic expression Eq. (j4"T|) . 

The formulas for a large T, Eqs. ff26|) and (JIT]) , used usually for computing 
the rate of the transition [191 I2Q] can not be applied to find the finite-size 
correction, because Taylor expansion of g(u)) in 

J duj g {u) ( 25in[ ^ T/2] ) 2 = 2*^(0) + 0(1/T), (49) 

leads diverging integrals over u to the non-leading terms. Hence we estimate 
the integral over the time directly. 

Each state in initial or final states of S[T] is defined by its position in 
addition to the momentum in field theory. The invariance under the trans- 
lation is broken slightly, and S[T] does not commute with the momentum 
operator 

[S[T],P]^0. (50) 

From Eqs. (T4*0l and (|50l) . both of the energy and momentum are not conserved 
in S[T]. Finite-size corrections of the transition rate and others are computed 
with S[T] of satisfying the boundary conditions and with the T-dependent 
expression Eq. fj4"5]) . We will see that particle states at ultraviolet energy 
regions in a relativistic invariant system, which couple with the S[T] in a 
universal manner, give the finite corrections. 
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2.2 Internal symmetry 

When a conserved charge Q commutes with the total Hamiltonian 

[Q,H] = 0, (51) 

and with a free part, H , 

[Q,H o ] = 0, (52) 

then Q commutes with the S-matrix at the finite-time interval, S(T), 

[Q,S(T)] = 0. (53) 

Hence the charge connected with an internal symmetry is conserved in S(T). 
Thus a state \ip) and the transformed state S'[T]|'0) have a same charge 

w> = (54) 

QS[T]\^) = qS[TW. (55) 

Hence the charge is conserved in the transition process of the finite-time 
interval also. 

If Q is a charge of non-compact group, its eigenvalue 

Q\q) = q\q), (56) 

is continuous and the eigenstates are normalized with Dirac delta function, 

(qi\q 2 ) = 2Ti5(q 1 - q 2 ). (57) 

Then the matrix element is written in the diagonal form in q, 

(q2\S[T]\q 1 } = 2rr5(q 1 -q 2 )S[T](q 1 ), (58) 

and the matrix element between any states is written with the reduced matrix 
S[T] 

[ dq 2 dq 1 F(q 2 )(q 2 \S[T]\q 1 )G(q 1 ) = 2tt [ dq 1 F(q 1 )S[T](q 1 )G(q 1 ). (59) 
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2.3 Unitarity 

The S[T] is defined with M0ller operators, Eq. fl38|) . and satisfies a unitarity 
relation, 

S t [T]S'[T] = 5[T]S f [T] = 1, (60) 
and an optical theorem 

i(T[T] - T[T] f ) = T[T]T[T] t , (61) 
S[T] = l + iT[T\. (62) 

The probability is preserved in S[T] and the imaginary part of the amplitude 
at T is determined by the total probability measured at T. 

S[T] is decomposed into the energy conserving term Ti[T] and non- 
conserving term T^T] 

S[T] = l + i(T 1 [T] + T 2 [T]), (63) 
then the unitarity Eq. (16 Op is written in the form, 

(l + i(T![T] +T 2 [T]))(1 -i(T}[T] +T 2 f [T])) = 1. (64) 
Hence we have 

-i^T} - Tl[T\) - z(T 2 [T] - T 2 f [T]) = 7i[T]T^[T] + T 2 [T]T 2 f [T] 

+ 7\ [T]T 2 [T] + T 2 [T]Tf [T] , (65) 

and 

- z(Tx[T] - T/fT]) = 7i[T]Ti [T], (66) 

- z(T 2 [T] - T 2 [T]) = T 2 [T]T 2 f [T] + Ti[T]T 2 [T] + T 2 [T]T/[T]. (67) 
The total transition probability from a state a is 

P= Yl l(/ 3 l T il«)| 2 + E l^l^la)! 2 , (68) 

where the energy conserving term is expressed by 

(0|Ti|a) = (/3|f 1 |a)27r5(^ - S a ). (69) 

It is noted that the unitarity connects physical quantities measured at T. 
Optical theorem proves that the imaginary part of forward amplitudes at T 
is written by total probability at T. Hence the life time at T, depends on T if 
the finite-size correction to the total probability is finite. The unitarity does 
not connect the probability at T with those at a different T. 
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2.4 Fluctuations 



Fluctuations of elementary particles are specified by Compton wave length 
Lcu> = — , where m is a particle's mass. L cw for various particles are 



Leu; — 



'2 x 10" 


7 


m] 


neutrino im v < leV/c 2 


4 x 10" 


-13 


[m] 


electron] 


< 1 x 10" 


-15 


[m] 


muon; 


2 x 10" 


-15 


[m] 


pi on; 


2 x 10- 

V 


-12 


[m] 


proton. 



(70) 



They are short range. 

In a relativistic field theory, a fluctuation of a field is expressed with a 
Green's function. A wave produced at a space-time position x\ and observed 
at another position x 2 is a Green's function, A(xi — x%) and is a function 
of x\ — £2 owing to a translational invariance. A(xi — £2) is written in the 
Fourier transformation 

^-^-^/^^"-'f^?- (71) 

where an integration in p° is taken along a complex path which avoids a 
pole of the integrand. From the residue of the pole, we have a component of 
on-mass shell waves of positive energy 

Ao(*i - x 2 ) = ^3 J i^/ Hxi ~ X2 \ E(p) = VP + ^ 2 - (72) 

Ao(xi — £2) is composed of (1) a singular part that is proportional to the 
Dirac delta function of the form S(X), A = {x\ — X2) 2 and (2) Bessel functions. 
The former one of (1) is narrow in A and long-range in \t\ — or \x\ — x 2 | and 
is called light-cone singularity. The light-cone singularity is characteristic of 
relativistic invariant systems. Since a one-particle energy of a mass m and 
a momentum p is E(p) = ^/p 2 + m 2 and approaches E(p) — > \p\ in a large 
momentum region, the phase in Eq. (1721) . p ■ {x\ — x?), cancels at a light- 
cone, |ti — £2! = l^i — a?2 1 i n the direction p in the infinite momentum. The 
wave becomes singular and real along the light cone and gives important 
contributions to a probability at a finite distance. The latter ones (2) are 
short-range functions which are specified by a length similar to Eq. ( !70|) and 
describe short-range fluctuations. 



16 



2.5 Boundary conditions 

Boundary conditions are necessary to determine a solution of a wave equation 
uniquely but is not explicitly written sometimes. In scattering processes, 
asymptotic boundary conditions are important [H]. For a scattering of a 
scalar field from an initial state \a) to a final state \/3), the states \oc) at 
t = — oo are constructed with free waves and the states \0) at t — oo are 
constructed also with free waves and satisfy asymptotic boundary conditions, 

lim (aW{t)\P) = (a\4 n \P), (73) 

t— >— oo 

hm(a|^(t)|/3) = (a|0Ll/3)- (74) 

t— >oo 

In the above equations, the field operators, <j>f{x), 4>{ n (x) , and 4>l ut (x) are 
expanded with a complete set of normalized functions f(x,t), and the coef- 
ficient <p*{x) is given in the form, 

tf{t)=i [ d 3 xr(x,t)d <P(x,t), (75) 



and the same coefficients 4>{ n (x) and (f>l u t( x ) are defined in the equivalent 
manner. 

To study a scattering at a finite-time interval T, asymptotic conditions 
are modified to 

hm (a\<pf(t)\P) = (a\<p{ n \P), (76) 

t— >~ T/2 

hm (a|^(t)|/3) = (a|0Ll/3)- (77) 

t^-+T/2 

In the above equations, the fields 4>{ n {x) and 4>lmt{ x ) an d the function f(x,t) 
satisfy the free wave equation and the states \a) and are defined with 
wave packets. Since the wave packets have finite spatial sizes and decrease 
fast at large \x— xq\ around a center xo, they ensure the asymptotic conditions 
at a finite T. It is important to satisfy the conditions even at a finite T to 
compute the finite-size corrections that are compared with experiments. We 
will see that in fact the finite-time interval corrections are obtained with wave 
packets. 
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3 Coherence length and asymptotic behav- 
ior of pion 

A charged pion is produced in a hadron reaction and decays to a neutrino and 
a charged lepton, after it propagates a certain distance. A pion is described 
with a wave field and its wave nature in the reactions is determined by the 
wave equation and the transition amplitude. We study hadron reactions 
using a Lagrangian of the pion and nucleon fields and analyze a position- 
dependent probability of detecting a pion at a finite distance. It is found 
that the position dependence of the probability is determined by the pion's 
mass and energy. 

Medium effects and detector effects are studied in Appendix and the 
pion's decay is studied in II. 

3.1 Pion emitted from a nucleon 

Pions and nucleons are an iso-triplet and an iso-doublet, and are expressed 

— * 

with fields <f)(x) and N(x). This physical system is described in term of the 
renormalizable Lagrangian, 

C = iV( 7 ■ p - m N )N + gN lb r ■ ${x)N + X - (dj^ - ^m 2 J 2 (x), (78) 

where and m n are masses of the nucleons and pions. A mass difference 
between a proton and a neutron and that of neutral and charged pions are 
ignored and SU (2) symmetry is assumed in most places. Second term in the 
right-hand side shows an interaction between a nucleon and a pion. Due to 
this interaction, a nucleon emits and absorbs a pion in intermediate states. 
These physical processes are treated by a renormalization prescription in a 
nucleon self-energy and others. The self-energy has no imaginary part at a 
mass-shell momentum since one nucleon does not decay and is stable. These 
particles in asymptotic states, which are defined at t — ±oo, do not decay 
due to the energy-momentum conservation. A nucleon at t — — oo does not 
transform into a nucleon and a pion at t — +oo. Nevertheless, one pion and 
one nucleon can exist for a short period as a virtual state. So this pion could 
be detected with a finite probability in a small space-time area using a small 
apparatus. Hereafter a space-time area where one pion is observed with a 
finite probability is estimated. 
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Let a nucleoli of a momentum p^. set at time t = U make a transition 
to a pion of average values of a three dimensional momentum p n at a four 
dimensional position (T W ,X W ) and other particles. A pion wave function of 
having these properties is w(p, X) of Eq. (j4]). This set of functions is complete 
set [21] and is equivalent to that of plane waves mathematically. Now the 
amplitude expressed with w(p, X) satisfies the boundary condition at T, and 
is a part of S[T], hence is different from the amplitude expressed with the 
plane waves that satisfies the boundary condition at T = oo. The probability 
of measuring the momentum and position with a small detector composed of 
nucleus of the transition process is computed with this amplitude expressed 
with wave packets. The wave packet size is a size of wave function that 
the pion interacts with. 

In a transition amplitude of a nucleon of the momentum p^ i at t — Tjv 4 
to a nucleon of the momentum p^ f and a pion, 

T = J d A x (Nf^ionlH^lNi),^ = gN l5 r-$(x)N, (79) 

the time and space are integrated over the region T^ < x° <T n , X 3 N . < x J < 
X^. The time interval is finite and this amplitude receives contributions 
from both of conserving and non-conserving states of the kinetic energy, 
from Eqs. (123]) and (HUj) . The former is a normal term and agrees with that 
computed with S-matrix of plane waves and the latter is a correction term 
that vanishes at T — > oo and is not computable with S-matrix of plane waves 
but with S[T]. 

The initial and final states are described by plane waves or the wave 
packet of the central values of the momentum and coordinate, and the width 
in the form 

\Ni) = \p Ni ,T Ni ), \N f ,pion) = \N f ,p Nf ;pion,p p i on ,X pion ,T pion ). (80) 
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The matrix elements of these states are defined in the forms, 

3 



l l 

77i jv \ 2 / rriN \ 2 1 



(27r) 3 \E{p Nf )J \E{p Ni )J ^Vi 
x uipN^uipNy^^-^^ (82) 



where 

3 - i 



A^- = (— V , p n = ( — ) , Vj : volume for initial state. 

V 7r / \2.£/ 7r (27r)' J / / 

In this paper, the spinor's normalization is 

U (P> s Mp, s) = 7 2^ m - ( 83 ) 

s 

In the above equation the pion life time is ignored. a n in Eq. ([8*1]) is the 
size of the pion wave packet. Minimum wave packets are used in most of 
the present paper but non- minimum wave packets give the equivalent result, 
which are presented in II. 

Substituting Eq. (IHTj) . we have the amplitude, Eq. (T~~~~) . in the form, 

N J dtdxdk e - *^* - ^/ - **)'"^ 5 ^/) -1 "^ 

= N J dtd£ n (2n) 3 S L @N t -PN f ~ 4) e < E (^/)- B (^)+ B ^)*+ i£; ^) T ^ 

x tZ(p JV/ )7 5 u(pi Vi )e-^ t ^-^ 3 e- < ^) Tir -^-^, (84) 

i 

m N rriN 1 ^ 2 



N = N o 

kP * \(2nrE(p Nf ) (2ir)*E(? Ni ) VV>, 

where S-^(x) is an approximate delta function of a finite size L, where L = 
\X W — AjvJ- The real part of the exponent becomes stationary at the center 
of Gaussian function, 

K = P-k- (85) 
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Combing this momentum with 



PN f = Pn, ~ K, (86) 

the f>N f satisfies 

PN f ~ PNi ~ Pit, (87) 

and the integration around this momentum has a finite contribution to the 
probability. Since momentum is conserved within uncertainty -^=, that is 
exactly conserved in the limit, — > oo. So the probability from this kine- 
matical region is regarded as an energy- conserving term that is computed in 
the standard method of plane waves and asymptotic boundary conditions. 
This probability is independent of a w and agrees with the value obtained 
by the ordinary S-matrix. Thus a small violation of the energy-momentum 
conservation does not modify the total probability of the normal component. 

Now there is another stationary momentum that is obtained from the 
imaginary part of the exponent. At large — and \X n — -XjvJ, the 
imaginary part of the exponent becomes large and the amplitude oscillates 
rapidly unless it satisfies the stationary condition, 



4-^ = 0, £ = {(E(p Nf ) - E(p Ni ) + E{k))t - E(k)T n + k ■ 

Ok-n 

The solution is, 



(88) 



v Nf t - v w t + v n T^ — Xn — 0, 

vn, = i^r, V* = ZZ l r Z'\ - (89) 



PN f Pn, - PNf 

E{p N fY n E(p N . -p Nf ] 



The second stationary momentum Eq. (189j) is different from Eq. ( 1861) and 
causes a new contribution that becomes important at larger T„- — T^. and 
iX-x — Xtf.\. The second stationary momentum is determined with the space- 
time position and k w could deviate from p n . Hence the violation of the energy 
and momentum conservation becomes large. Thus two stationary momenta 
have different properties. First one does not depend upon the time but the 
second one does and they lead different phases. The violation of the energy 
and momentum conservation is small in the first one but is rather large in the 
second one, hence the first one corresponds to the normal term of Eq. (HT|) 
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and the second one corresponds to the finite-size correction 5f of Eq. ( 141 p . 
The two terms of different phases lead the probability to have a large finite- 
size correction. We will see next that the finite-size correction is computed 
rigorously and has a universal property. 

3.2 Expressing probability with correlation function 
and light-cone singularity 

The total probability is expressed in the form 

1 1 



P = dX. 



where P is written from Eq. ( 179 p . 



P= dxidx 2 A(xi,X2)e ^ 2 



J2i ^{Xi-Xir-viU-Tir)) 2 iE n (pi r ){tl-t2)-ip 7 r-(x 1 -X2) 



(90) 



(91) 



with a correlation function 

A( Xl ,x 2 ) = J Ndp N ^p Ns J Ni )e-^- pN f Hx '- X2 \ (92) 

d(p Nf ,PNi) = \ ^(^PN f h5u{m)Tu{pN f h5u{p Nl ), (93) 



spin 



N 



m 



N 



(2n) 3 E_ 



Taking a sum of the final spins and an average over the initial spin, we have 
d(pN,,PN f ) = (PN t ■ PN f + m 2 N )/2m 2 N , (94) 
and the correlation function 
1 



A(x 1 ,x 2 ) 



(2k) 



dpNf ,2 i \ -i(pN -PN f )-(xi-x 2 ) 

'm N +p Ni -p N )e va, i 



E 



N f 



d 4 p9(p°)lm 



p 2 — mjj + ie 



(% + PNt 
1 

p 2 — m 2 N + ie 



p^ e ~ i iPN i -p)-(x 1 -X 2 ) 

(m 2 N +p Nr p)e ip < x '- X2 \ 
(95) 
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The integral over X n in Eq. ( 19 Op gives 



3.2.1 light-cone singularity 1 

Here the formula for a relativistic field 



1. 



(96) 



—5(X)e(5t) + f short 

47T 



J d 4 q6(q°)5(q 2 + m 2 )e iq - Sx = (2tt) 3 z 

9(-X) jiVi (mv^) - ie(<Jt) Ji (mv 73 !) } 



im 



- 0(A) ^^ ifi f tWa) , A = <fe 2 , <Jt = 
47r 2 VA v / 



(97) 



where iV x , J 1( and i^i are Bessel functions, is substituted into Eq. ( 195|) . Then 
A(xi,X2) is written as a sum of a singular term and other terms expressed 
by Bessel functions of A = (x\ — x 2 ) 2 



A(.r,..r,) = ,-iPN^i-x 2 ) j ^^(h - t 2 ) + Bessel functions^ 



-i<t> N (St) _ 



47T 



5(X)e(8t) + Bessel functions, 



(98) 



>N.. 



, = (E Ni - \p Nt \)5t, 5t = t 1 -t 2 . 



The first term in the right hand side of Eq. ( 1971) is the most singular term 
and the second and third terms have singularity of the form 1/A around 

A = and decrease as e _m V / ^' or oscillates as e %m ^^ . The singular functions 
and regular functions behave differently and are expressed in Fig. [2] for one 
space dimension. The singular functions are non- vanishing in narrow regions 
around the light cone and the regular functions have finite values in a small 
area around the origin. Since the light cone is extended in macroscopic area, 
the light-cone singularity makes the correlation function long-range. Thus the 
correlation function 1\ becomes long-range only along the light-cone region 
and decreases exponentially or oscillates rapidly in other directions. In the 
above equation, "Bessel functions" are oscillating or decreasing fast with A 
and this property is sufficient to know a large T behavior of the probability. 
In latter sections, concrete forms of these functions are obtained and their 
properties are studied. 
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Fig. 2: The space time region of the correlation function is shown. The region 
I around the origin corresponds to short-range correlation where t ~ x ~ 0, 
while the region II corresponds to long-range correlation where t 2 — x 2 ~ 
and both t and x can be macroscopic. 

Eq. ( 198]) is substituted into Eq. (l9Tj) and after a tedious calculation which 
will be presented in II, the probability is written in the form, 

P = I " dhdUa^ ^^ef^e^W^iW + 5esseZ functions, 
Jt Nx 2 \5t\ 

4>n = ujt, u w = (E v — IpVI), ^ = h - t 2 , (99) 

where we use the notation t = x° for the sake of simplicity. The probability 
derived from the first term of Eq. (19"9"1) is proportional to the following function 
of cuT expressed as 

T e iu,6t x / . 1 - cos x 



Tg{T,u) = i I dt 1 dt 2 -—e(5t) = -2 (Sin x ) , (100) 



o 



x = ujT, Sin x = I dt . 
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which satisfies 

JLg(T,u)\ T=0 = -OJ, (101) 

g(oo,u u ) = —n. (102) 
Subtracting the asymptotic value, we define g(T,u u ) 

0(T,w) = g(T,u) +g(oo,u), (103) 

which satisfies 



d . / sin(x/2) 



2 



_$(*) = -4 , (104) 

ax \ x 



and oscillates rapidly at large x. Averages become 



g(x) = -. (105) 
x 

Thus g(T,uj) behaves as 

g(T,u)n^, T = -, (106) 
i to 

and is given in Fig.[3j 

The probability derived from the second term of Eq. (1991) is a constant. 
Thus the probability is expressed in the form 

P = C' g(co,T) + P , 

co = (E n - \p„\) - (E N . - \p Ni \), T = T w - T Ni . (107) 

g(uj, T) is a function of wT and is inversely proportional to uT at a large loT. 
Here to is determined by the pion's mass and energy and the nucleon's mass 
and energy and is not very small. Hence, g(u)T) vanishes at a macroscopic 
T. The integrand of the probability in Eq. (1991) oscillates rapidly in St with 
a non-small angular velocity u, and the probability becomes constant fast. 



3.2.2 light-cone singularity 2 

We should note that the integrand of Eq. (195]) around a momentum region of 
(PN f ~ PNi) 2 — does not oscillate at A = and becomes real. Those waves 
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Fig. 3: The function g(T, oS) are given for various particles. The mass m v = 1 
[eV] is used for the neutrino in this Figure. Since the electron and muon are 
massive, g(T,u) are negligibly small. The values for the pion and other 
hadrons are smaller than that of muon. 



that are added constructively form the light-cone singularity of the form 
that does not accompany the oscillating function e~ l ^ N ^ St \ This function 
is extracted in the following way. We change the momentum from p to 
q — PN f — Pi and have the expression with the new variable q 



A(x 1 ,x 2 ) 



d 4 qe(p% - q°)Im 



1 



{2m 2 N -p Ni - q )e~^-^\ 
(108) 
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Next the denominator is expanded in the form, 
1 1 



(<? - PnJ 2 ~ m 2 Nf + ie q 2 + m 2 N . - 2q ■ p N . - m 2 Nf + ie 

1 

q 2 - 2q ■ p N . + 5m 2 N + ie 
i 



J2[ 2q - p ^ dSm %) q 2 + Sm 2 N + ie 



(109) 



where a small difference between m^. and mjsr, which has been ignored so 
far is included here and 5m 2 N = m 2 Ni — m 2 N is the mass-squared difference 
between a proton and a neutron. We substitute this expression to A(xi,x 2 ). 
From the integration of q° < in Eq. fllOSp . the light-cone singular term 

A(xi,x 2 ) = —6(X)e(St) + Bessel functions, (HO) 

is derived from the lowest / = term of Eq. (I109p . Eq. (IllOp is convenient to 
compute the large St behavior of A(xi,X2). A(xi,X2) has also non-singular 
terms. This method of extracting the singularity is valid when the infinite 
series converges. The integrations of the series in the expansion of the de- 
nominator converge in the kinematical region, 

2Pn z ■ Pn < 5m N 2 . (Ill) 

Hence the present method is valid and the light-cone singularity exists in this 
narrow kinematical region. The convergence of the series will be studied in 
II in detail. 

Eq. (IllOp is substituted into Eq. (I9~T|) and after tedious calculations, we 
have the probability in the form, 

* d^M^^rreiSty^W + others, 
2 \ 5t \ 

= CoK)£K,T) + P . (112) 

In Eq. (I112p . P is the normal term and the first term is the finite-size cor- 
rection. The correction is the product of the universal function g(u n , T) that 
is independent of the wave packet and Co^a^). 
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Eq. ( |112p has almost the same form as Eq. (I107p but a different angular 
velocity u n is used. co n is smaller than u of Eq. (11071) . Hence this is more 
convenient to study the large T behavior than the former one since it has the 
slowest oscillation. At large \St\ region, the frequency in time is given from 



In the last calculation, the large momentum expansion E{p) = \p\ + S^r was 

2 

made. The pion's mass m n , 139 [MeV/c 2 ], makes the angular velocity 
small only in an extreme high-energy region. Thus a coherence length of a 
pion emitted from a nucleon is Iq = ^ ZL , which is microscopic due to the 



large mass, and the pion in the kinematical region, Eq. fillip can be observed 
inside this length. If the pion's mass m n were 1 [eV/c 2 ], then l would have 
been macroscopic. 

3.3 Pion from NN collisions 

A pion produced in a nucleon collision shows the same behavior as a pion from 
an isolated nucleon of the previous subsection. Since a mass-shell nucleon 
lives the longest period of time compared to a virtual nucleon state, a pion 
emitted from a real state of nucleon has the longest correlation. This pion 
is expressed with external lines in Feynman diagram Fig. 4 . We study this 
process hereafter in this section. Using a scattering amplitude of a proton 
with a target, T(px f , ■ ■ ■ ,p n ;Xr), an amplitude of detecting one pion of a 
momentum p n at a position is 



(Pxdknwipn - K; X^)e l(pN f +kn PN ^' x u(p N )^u(p Ni )T(p Nt} ■■■ ,p n ; X T ), 



where p^. and pi,P2, • • • are momenta of the initial and final state, and X? 
is the position of the target, and x\ is the coordinate where the pion in the 
external line attaches to the nucleon. We study a process expressed with 
a diagram of Fig. 4. Using a scattering amplitude, a correlation function is 
expressed in the form, 



<f>{St) = {E v -\^\)6t 




(113) 



:2 



(114) 



C^X^tPk) = J dp Ni P(X- K ,p n ; 
w{p- k, X) = N' e~ a $-v) 2 - ik - 




(115) 
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Fig. 4: A Feynman diagram of a pion production in NN scattering, where the 
pion has the longest correlation length. The pion emitted from the external 
nucleon has the longest correlation and is computed in the text. 

where P(X n ,p n ;pNi) is the probability of a pion from a nucleon of the mo- 
mentum given in the previous subsection. Thus the behavior ofC(X n , p n ) 
is equivalent to that of P(X n ,p n ;pN i ) and has a length l of Eq. (fT3l) . The 
length Iq = is much longer than the de Broglie wave length, A = A, but 
is microscopic value m ordinary high energy experiments. 

Thus the pion's coherence length Iq is of microscopic size and experiments 
of detecting pions have no finite-size corrections. Pions are produced in 
small area at hadron collisions with certain probability and in the region 
IX^ — Xt\ ^> ^0) they behave like particles. 
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4 Summary 



The probability to detect a pion emitted from a nucleon is given in Eq. (11121) 
and is decomposed to the normal term Pq and the finite-size correction 
Co((j)g(ujT). The finite-size correction is expressed by a universal function 
g(cuT) in relativistic systems, and is stable under the change of parameters. 
The magnitude Co(er) depends on a and increases with a. The finite-size 
correction is small for a pion and large for a light particle. From the position 
dependence of the correction, new physical quantities may be lead in light 
particles. 

Thus scattering and decay processes of finite-time intervals are different 
from those of the infinite-time interval and have various unusual properties. 
A wave function started from an eigenstate of Hq and evolves with an inter- 
action Hamiltonian H\ is a superposition of states of continuous spectrum 
of H that includes the state of infinite momentum at a finite time. Despite 
the fact that the total energy defined by H is always conserved, the kinetic 
energy defined by the free part H Q is not conserved at a finite-time interval 
due to a many-body interaction effect. 

At the infinite time t — > oo, only the states of Ep = Eq couple to the wave 
function. Consequently the asymptotic values are determined by the states of 
conserving kinetic energy. At a finite time, on the other hand, various states 
of different kinetic energies couple. These states of continuous kinetic energy 
give a non-uniform behavior to many-body wave function and a diffraction 
phenomenon. 

The amplitude and probability to detect particles at a finite distance re- 
flect the unusual behavior of the wave function at a finite time. They satisfy 
boundary condition at a finite-time interval T, and are formulated with the 
S[T), S-matrix at a finite T, expressed by wave packets. Wave packets are 
known for a long time to be a powerful tool for formulating scattering am- 
plitudes in rigorous manners. But they were replaced with plane waves in 
practical calculations, because the finite-size corrections have been thought 
ignorable in high energy regions so far. We show in this paper that the wave 
packets are not only useful for formulating the scattering amplitude but also 
can be used for computing finite-size corrections of the probability to detect 
the particles in quantum mechanical processes. Particle's energy is given in 
a relativistic system as Ep = a/p 2 + m 2 and approaches \p\ at the ultra- 
violet region. Hence the corresponding waves have almost the light velocity 
and cause the light-cone singularity to correlation functions. The light-cone 
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singularity helps to compute the finite size corrections at long range region. 

The finite-size corrections to the probability to detect the pion produced 
in hadron reactions are studied as an example in I. We found that the coher- 
ence length for the wave phenomena to appear is small for the pion, and the 
correction is negligible in ordinary experiments. 
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Appendix A Wave packet shape 

Scattering matrix is expressed by wave packets in rigorous manner [2] . Wave 
packets are important, therefore, but they can be replaced with plane waves 
in practical situations in the asymptotic regions where the final state can 
be treated as free particles. In non-asymptotic regions where the final state 
retains wave nature, the wave packets are necessary to compute the position 
dependent quantities. The wave packets may be unfamiliar to some readers 
and we give a summary of wave packets. 
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A particle of a finite coherence length is described with a wave packet, 
which is extended in momentum and localized in position around the centers 
(p, X) and satisfies a corresponding wave equation. Although its precise 
shape may be unknown generally in real experiments, the physical quantities 
of universal properties that are independent from the details of wave packets 
are important. 

First, the wave packets are localized in the momentum and position 
around their centers [23j I2U [251 [26] . The whole set of the wave packets 
becomes complete set [2Tj . 



independent of the shape parameter (3 which is discussed later. 

Second, we study the wave packets that preserve the discrete symmetries 
such as invariances under space and time inversions, which are quite natural, 
since the origin of the wave packet is the interactions of the particle with 
matters in detectors. The wave packet sizes are estimated in the next section 
under this consideration. We study, hence, the wave packets which are su- 
perpositions of the plane waves around the central momentum with a weight 
function that has the same property under these transformations, 



where the momentum p is the central value of the momentum. In the present 
work, k is used for the integration variable and p is used for the central value 
of momentum. 

Since the time reversal invariance is satisfied in the quantum electrody- 
namics and weak interactions in the lowest order, the wave packet which 
is invariant under the time inversion is most important. Under the time 
inversion, the coordinate and momentum variables are transformed into 







x — > x, k —>• —k, p — p, t —t, 




and the plane wave is transformed to its complex conjugate, 



i(Et—k-x) . —i(Et-k-x) / i(Et—k-x)\* 



(119) 



So when the weight satisfies 




(120) 
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the state described by the wave packet is the superposition of the waves 
which are transformed equivalently under the time inversion, 

dkw(k;p)e i{Et -^ -> J dkw{-k--p){e i{Et - t£) )* = (J dkw(k-p)e l{Et ~ l 

(121) 

Next we study the space inversion. Although neutrino violates this sym- 
metry, the wave packet may preserve. Under the space inversion, the coor- 
dinate and momentum variables are changed into 

x — > —x, k — > —k, p — > — p, (122) 

and the plane wave is changed into 

e i(Et-k-x) e i(Et-k-x) ^ (123) 

So when the weight satisfies 

w(k;p)^-w(—k;—p)=w(k;p), (124) 
the state described by the wave packet is transformed in the following way, 

J dkw{k]p)e i{Et ~ tiZ) J dkw{-k;-p)e i{Et - ts) = J dkw(k;p)e i(Et - ls) , 

(125) 

in the same way as the plane wave. 

The invariance under the time reversal, Eq. f)120p . is a strong condition 
and leads the important result. We study wave packets that satisfy Eq. (I120p . 
The simplest form of satisfying this property is the Gaussian wave packet 
\p,X,Po) 

|£X,/3 ) = f d k e -^-v)\i^-T)-H^-x))^ (126) 

(2tt) 



) 2 



where the parameter a shows the size of the wave packet in the coordinate 
space and iV is the normalization factor. We write a wave packet of a pa- 
rameter (3 as \p, X, (3). The Gaussian wave packet is the state of (3 = O and 
non-Gaussian wave packets are the states of (3 ^ (3q. The calculations using 
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the Gaussian wave packet are presented in most places for the sake of sim- 
plicity, except the derivation of the large distance behavior of the neutrino 
wave function, which will be made in II. 

The normal physical quantity of microscopic physics is obtained from 
ordinary scattering which has no dependence upon the distance or time in- 
terval between the initial and final states. This is because the size of the 
microscopic system is so small that the experimental apparatus is regarded 
infinite and the boundary conditions of ordinary scatterings which are defined 
at t — ±00 of plane waves are suitable. Hence the amplitude and probability 
are defined using plane waves and the boundary conditions at t = ±00 ensure 
the independence of the probability of the distance and particle's coherence 
length. 

From Eq. f ll 16 j) . the total probability of the transition A — > B + C is 
independent of the parameter (3 of C 

]T \(C;p c Jc,(3;B\T\A;p A ,X A ,T A )\ 2 

PcSc 

= mp c ,Xc,f3o;B\T\A;p A ,X A ,T A )\ 2 , (127) 

vcSc 

which also agrees to that of the plane waves. The probability of the finite 
distance is defined by restricting the center positions X inside a finite spatial 
region V. The completeness condition for the function in this region V is 
then 

\p,X,f3 )(p,X,f3 \= \p,X,(3)(p,X,f3\, (128) 

P,X<V P,X<V 

and the probability satisfies 

£ I (C; pc, X c , (3; B\T\A; p A , X A , T A )\ 2 

pc,x c <v 

= mPc,X c ,f3 -B\T\A-p A ,X A ,T A )\ 2 , (129) 

p c ,X c <V 

and is also independent from the /3. We study the probability that is the 
average over a finite energy region V p 

£ \(C-p c ,X c ,(3;B\T\A-p A ,X A ,T A )\ 2 , (130) 

Pc<V p ,X c <V 
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and extract the physical quantity. For the genuine physical quantity, the 
value should be determined uniquely. For the neutrino experiments, the 
energy uncertainty is of the order of 10 percent of the total neutrino energy. 
This magnitude is the same order as that of the minimum uncertainty if the 
neutrino energy is 1 [GeV]. So the minimum wave packet is applied directly 
in this case. The probability for a larger energy uncertainty is computed 
using the probability of the small energy uncertainty. The dependence of 
the neutrino probability on the distance will be shown to have the universal 
behavior. 

The S-matrix at a finite-time interval T, S[T], depends on the bound- 
ary condition at T, which is determined by the wave packet size a of the 
final state. Hence the probability from S[T] and finite-size correction are 
independent of parameter but depends on a. It is found in fact that the 
finite-size correction has the factorized form of g(uiT) and C(a). Only the 
latter depends on a. 

Appendix B Wave packet size 

We estimate the wave packet size of a proton first and those of a pion next, 
following a method of our previous works [2TI [22] . 

B-I Proton mean free path 

A mean free path of a charged particle is determined by its scattering rate 
with atoms in matter caused by Coulomb interaction. An energy loss is 
also determined by the same cross section. Data on the energy loss are 
summarized well in particle data summary [7] and are used for the evaluation 
of the proton's mean free path. 

The proton's energy loss rate at a momentum, 1 [GeV/c], for several met- 
als such as Pb, Fe, and others are 



hence we have the mean free path of the 1 [GeV/ c] proton in the material of 
a density p, 



dE 

dx 



1-2 [MeVg _1 cm 2 ], 



(131) 



L 



E 



1 [GeV] 



50 - 100 [cm]. (132) 



proton 



x p (1 - 2) x 10 [MeV g^cm" 1 ] 



dx 
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At a lower energy of the order of 0.2 [GeV/c], the energy loss rate of the 
proton is about 10 [MeVg _1 cm 2 ] and the mean free path is 



£ P roton = 10 [cm]. (133) 

A wave which describes a proton maintains coherence in matter for a finite 
distance of the mean free path, hence this wave is approximately described 
by a wave packet that has a size of the mean free path. We use the mean 
free path for a wave packet size of the proton A/o" proton , 

V^proton = -^proton- (134) 

When a proton of this size is emitted into the vacuum or to a dilute gas 
from matter, the wave keeps the same size. The proton that is moving freely 
has a constant size in vacuum or dilute gas. The size varies when the proton 
is accelerated. If the potential energy V is added to the proton of momentum 
Pbefore, then the final value of the momentum becomes p a ftcr and satisfies 



before + ™* + V = yj + m\ (135) 

From Eq. (II 35ft variants of the momentum satisfy 

V before X ^Pbefore = ^after x Rafter, 

Pbefore Pafter /i oc\ 

^before = , ==, ^after = , ==• UoOj 

V Kefore + m V Rafter + m 

Hence the size of a particle, \/<7 before , which is proportional to the inverse 
of Sp he{ore , becomes A/o" after after the acceleration from a velocity v before to a 
velocity v after- The wave packet size is determined by the velocity ratio, 

V^aftcr = V^before X — (137) 
^before 

The velocity is bounded by the light velocity c, and a velocity ratio from 
1 [GeV/c] to 10 [GeV/c] is about 1.2 and that from 0.2 [GeV/c] to 10 [GeV/c] 
is about 5. Hence the proton of 10 [GeV/c] regardless of the energy in matter 
has the mean free path 

v^proton ~ 40 - 100 [cm], (138) 
in vacuum or a dilute gas. 
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B-II Pion wave packet 



Wave packet size of pions which are produced by a proton collision with target 
nucleus is determined by the proton's initial size Eq. dl 38[) and a target size 
10 -15 [m], which is negligibly small. A pion is produced while the proton wave 
packet passes through the small target by the strong interaction, hence this 
pion has a size in temporal direction of the proton wave packet. Hence the 
size of pion wave packet, A/o" pion , is given from that of the proton, y^proton' 
in the form 



^proton V ^pion / — ^pion / — / — /_. „„\ 

— ~ ~~ — ~ J V^pion ~~ ~ V°"proton ~ V ^"proton- I ■'■""J 

^proton ^pion ^proton 

In relativistic energy regions, particles have the light velocity. Consequently 
from Eq. f |138j) . pion's wave function of 1 [GeV/c] or larger momentum has 
the size 

v^ pion « 40 - 100 [cm]. (140) 
We use this value of Eq. f )140p as the size of the wave packet 

V 7 ^ = V^pion, (141) 

in II. 

In vacuum and dilute gas, pions of the size Eq. (j!41j) propagate freely. 
From Eqs. ( 11321) . ( j 138ft . fll40|) . the proton and pion have the sizes of the order 
of 50 - 100 [cm]. 



Appendix C Neutrino on target: neutrino wave 

packet 

A neutrino interacts in matter with a nucleon or an electron which are con- 
stituent particles in a bound atom and are expressed with wave functions of 
finite sizes. So a wave function expressing the final state of neutrino has a 
size of nucleus or atom, which is written with a wave packet. A size of wave 
packet for an observed neutrino, therefore, is not determined with a mean 
free path but with a size of the target in its detection process. They are 
either a size of a nucleons in a nucleus or that of an electron in an atom. Nu- 
cleus have sizes of the order of 10 -15 [m] and electron's wave functions have 
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sizes of the order of 10 11 [m]. So neutrino wave packet is either 10 11 [m] or 
[m]. 

Interactions of muon neutrinos in detectors are 



lO" 15 



+ e~ -> e~ + Up, 

+ e~~ -> H~ + "e, 
+ A ->• fi~ + {A + 1) + X, 
v^ + A^v^ + A + X, 



(142) 
(143) 
(144) 
(145) 



hence the size of the neutrino wave packet a/ov in processes ( 1142 j) and (11431) 
is of the order of 10 _11 ,[m] 



lO" 11 [ml 



(146) 



and the neutrino wave packet y/a\, in processes (j!44p and (I145j) is of the order 



of 10 



-15 



m 



-15 



m 



v^ = io- 

Interactions of electron neutrinos in detectors are 



v e + e — > e +f e , 
v e + A -> e~ + {A + 1) 
u e + A-te + A + X. 



(147) 



(148) 
(149) 
(150) 



The neutrino wave packet y/a^ in processes (11481) is of the order of 10 -11 [m], 
Eq. (I146p . and the neutrino wave packet ^[d\, in processes (1 149ft and ( 11501) is 
of the order of 10 -15 [m], Eq. (j!47p . They are treated in the same way as the 
neutrino from the pion decay. 

From Eqs. (1146 j) and (I147p . the neutrino have the wave packet sizes of the 
order of 10 -11 [m] or 10~ 15 [m]. 

For various nucleus the sizes are estimated from the nucleus size, ml 



as 



(5.2/ml;C 
6.35/ml;0 
14.3/m^; Fe 

{18.9 /ml; Pb. 



(151) 
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The probability to detect neutrinos is expressed with the wave packet for the 
detecting particle, which is determined by the size of nucleus in the detector. 
In this respect, the neutrino wave packet of the present work is different from 
some previous works of wave packets that are connected with flavor neutrino 
oscillations ^3I2S1I2SII3QIEI1E21I331E3], where one particle properties of 
neutrino at productions are studied. 

Appendix D Charged particles on targetrwave 

packets 

Charged particles are detected by their electromagnetic interactions with 
atoms in matters. The electromagnetic interactions are mediated by mass- 
less photons and the forces are long range and are much stronger than the 
weak interaction. Hence successive interactions with many atoms, which are 
correlated quantum mechanically each other in solid, give signals of the par- 
ticles. Thus the wave packet sizes of the charged particles would be much 
larger than the size of an atom but the size of many atoms. It would be 
reasonable to assume that the size is semi- microscopic, some number of the 
order of one times 10~ 10 [m]. This size might agree to those that have been 
considered before in textbooks [23J I2U E5J ES]- Although the wave packet 
sizes of the charged particles are much larger than those of the neutrinos, 
the diffraction components have the magnitudes, g(u c harge, T), which are ex- 
tremely small due to their large masses, and vanish at the macroscopic dis- 
tance. Consequently the finite-size corrections may be negligible for charged 
particles. 

D-.l Dependence on wave packet size 

It is known that the probability at T = oo does not depend on the wave 
packet size. The result of the present paper in fact shows that this probability 
is independent of the wave packet size. Now the finite-size correction is 
different and depends on the wave packet size. This peculiar property of 
the finite-size correction is one of the important results of the present paper. 
This correction has the universal properties, which were computed with S[T] 
defined by the wave packets. 
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